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ABSTRACT 

The American Association of State Highway and Transportation Officials (AASHTO) provide 5 methods for distributing highway 
superelevation (&) and side friction (f). AASHTO’s Method 5 uses traditional unsymmetrical curves and focuses on speed variations. 
The objectives of this paper are (1) to develop the superelevation distribution (@) that improves the design consistency of highways 
based on the safety margin, and (2) to find the best superelevation distribution for a specific highway using the two unsymmetrical 
vertical curves (EAU and SAU curve), and the optimization model. This paper applies AASHTO’s Method 5 algorithm and Easa’‘s 
equations on distribution of the superelevation (&) and side friction factor (f). The paper also describes © and f distributions based on 
the optimization model. The proposed methods improved highway safety by creating the geometric design consistency using the 
safety margin methodology. The safety margin is defined by the difference between the design speed and maximum limiting speed. 
The proposed methods improve highway design consistency based on safety margins and are superior to AASHTO’s methods with 
examples illustrated. 


Keywords: Highway design; Design consistency; superelevation distribution method; Optimization model. 
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1. INTRODUCTION 

Over half of the fatalities connected to collisions on rural highways happen on curved sections of the highway (Lamm et al., 1992). 
Therefore, the curved areas and their corresponding transition sections represent the most critical areas for safety improvements. 
One of the main methods for improving highway safety is to create geometric design consistency in the highway. Design 
consistency guarantees that the geometric elements of the road act in a coordinated way so that the performance required for safe 
navigation of the road is consistent with drivers’ expectations (Lunenfeld & Alexander, 1990; Krammes et al., 1995). 

Highway design consistency is evaluated through three main factors: the driver's performance, operating speed, and safety 
(Gibreel et al., 1999). Driver's performance is affected by his/her anticipation and workload, and also by road features such as 
interchange design and aesthetics. 

Naatanen and Summala (1976) proposed that drivers want to control the "safety margin," something that has been established 
for many areas within the engineering field. The safety margin of a road depends on its system of curves and may require 
differentiation between "the critical and real lateral acceleration, friction, or speed" (Summala, 1988). 

Nicholson (1998) however proposed a simple way of measuring the safety margin. Nicholson defined safety margin as the 
difference between the design speed and the speed limit. The greatest permitted design side friction is dependent on the driver 
comfort and the design speed. Nicholson (1998) applied this concept to a definite set of curves to evaluate the different AASHTO 
methods for superelevation distribution on highways. The author demonstrated that Method 1 (linear) is inferior to Method 5 
(curvilinear) based on the mean, standard deviation and coefficient of variation of the safety margin. Furthermore, the author 
emphasized the need to evaluate design consistency on the highway using the safety margin methodology. Using various data sets 
Nicholson (1998) and Easa (1999) have shown that AASHTO’s Method 5 is the best method. 

Easa (2003) developed an optimization model for the superelevation distribution that maximizes the design consistency of the 
highway based on the concept of safety margin. The quadratic model (aggregate analysis) findsthe best superelevation distribution 
for certain highway data using the entire design space of superelevation. It was evident through this implementation that the best 
method relies on the geometric characteristics of the highway. This was the motivation behind the research presented in this paper. 

Two unsymmetrical vertical curve equations have been developed by Easa (1994; 2007). The first one is the EAU (Equal Arc 
Unsymmetrical) vertical curve that contains two equal arcs with uneven horizontal projection for tangents (Easa, 1994). The second is 
the SAU (Single Arc Unsymmetrical) vertical curve which has a single arc, acquires the cubic function form, and lies above the flatter 
arc and under the sharper arc for the EAU vertical curve making it smoother (Easa, 2007). 

This paper provides two proposed methods and an optimization model for assessing the superelevation distribution that 
improves design consistency of highway, based on the safety margin. The proposed methods find the best superelevation 
distribution for a specified highway data using the two unsymmetrical vertical curves (EAU and SAU curve). 

The developed optimization model determines the best distribution of superelevation, where such an optimization model has an 
objective function subject to parameters. The optimization model does not only affect the difficulty of determining the best 
superelevation values utilizing trial and error, but additionally, it guarantees the most advantageous solution. Consequently, the 
model scans all probable values of superelevation within the complete area between Method 2 and Method 3 (Fig. 1) to find the 
best superelevation rate for specific highway data. 


2. AASHTO METHODS 

Whenever a vehicle moves on a curve within the design speed, the centrifugal force is counteracted by the superelevation with the 
side friction between the tires and pavement. Through the fundamental laws of mechanics the primary relationship of function ona 
curve can be written as in Eq. 1 (AASHTO 2004). 


etf=— (1) 


where e#=superelevation (in percent); f=side friction; ¥=design speed of the curve (or the average running speed of the curve) (m/s); 
g=force of gravity (9.81 m/s*); and R= radius of the curve (m). ¥7/(gR) is called the centripetal ratio. According to AASHTO (2004) 
when traveling at the precise design speed, the minimum turning radius (Ry) can be resolved such that both theé.caNd foooe are 
chosen: 


GoorET v= 


Poin = (0001 Cyryt finns) (2) 


where, ¥ is the design speed of the curve (km/h) and@yncx ANd fingyare the maximum permitted design of superelevation and side 


friction factor, respectively. 
For radii larger than Bynn, the design ofe and f are below their maximum rates. For a known design speed (@mex, fimax, along 


with Ryn), AASHTO provides five methods to counteract the centrifugal force through the utilization of ¢ and f (AASHTO 2004). 
These methods are described as follows (see Fig.1): 
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In Method 1the superelevation and side friction factor are directly 


g related to the inverse of radius R (otherwise the curve degree). That 
2 means estarts from zero and changes linearly. At R= = (curve 
(2) 2 degree equals Zero) into @.5. at R=Ryi, (Maximum curve degree), 
3 f also change linearly much the same way as ¢ does. Nicholson 
§ (1998) mathematically expressed Method 1 as: 
Bonus 7 
Oy = emax’ fe “5a g omow 
(Rin < R < =) (3) 
: where the subscript of # and f represents the method. Method 1 
2 can also be expressed as: 
(bo) § 
te a. - Cmax e - = Fax e : 
3 1 (Cena finad GR! !? (Bonar Fina GR 
ah (Rmin < R < @) (4) 
0 U/Ranin ‘ E ne ; 
UR In Method 2when a vehicle is driving at the design speed the 
- ae ee trif | f ted i irect ti to th 
Figure 1 AASHTO methods and the models of distributing pr alder oe ca peaanaebeetcge Pe counteracted e eine propane enveie 
, , aaa ‘ side friction upon curves of those needing fixe» For the sharper 
superelevation: (a) superelevation distribution and (b) side- oe a ; 
fictiondietibution curves, f stays at fing» while gincreases in direct proportion to the 


increase in curvature until ¢=€..... Nicholson (1998) showed that 
Method 2 can be mathematically expressed as: 


a> = - Fax’ f2= Fox (min < PR < Ry) (5) 


€:= 0; f2=G) Ryo < R < @) (6) 


The radius least relying on side friction is: 


Rfo = 55 — 7) 


In Method 3 for a vehicle driving at the design speed, the centrifugal force can be counteracted in direct proportion to the 
superelevation for curves that do not reach @.». For sharper curves,@ stays at @n.. while fchanges in direct proportion to the 
increase in curvature until f=fp,.¢». Nicholson (1998) mathematically expressed Method 3 as: 
ve 

@3= max fa= ge | max (Rmin < R < Rey) (8) 
And 

p= 
83 = ga’ fz = 9; (Res <R<« oo) (9) 


The radius that is the smallest relying on superelevation is: 


v= 


Rep = 4am (10) 


In Method 4, the average running speed is utilized as opposed to the design speed. Furthermore, the superelevation and side 
friction factor differ with the radius in an inverse curvilinear manner, providing values between Methods 1 and 3. 

Method 2 is modified using four methods to determine the superelevation and side friction. The second one only involves zero 
superelevation on infinite radius curves. In addition, the second method can also be modified to guarantee that a minimum 
superelevation is expressed. Nicholson (1998) shows how the modified method 2 can be explained mathematically as: 


v= 
82m = oR - Fax’ fom = Free ' 


(Bmnin <P < Bang) (11) 
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v= 
Crm = Smin/ fm = gu Smin 


(Rime < R < ©) (12) 


where 2m = "modified method 2." The smallest potential radius is used and the minimum superelevation applied is @y.23' 


v= 
me ~ F (Gm Fnae) 


R (13) 


In Method 5, the side friction and superelevation differ in an inverse curvilinear fashion with the radius, providing values between 
Methods 1 and 3. Easa (1999) shows how Method 5 can be explained mathematically as: 


ve 
Rp = 9 Smas (14) 


Machen 
fa = Femme WR < 1/ By (15) 
&. Ap;-# & a—&, 
fs= Frmase| a aa dd Saee (—ae ) 1; V/R>1/Rpy (16) 


ve 
as = To fs (17) 


The AASHTO guide has subjectively assessed these methods. The linear distribution for Method 1 is perfect when all vehicles are 
running at a constant speed and regardless of whether the movement of the vehicle is upon the tangent, the flat curve, or a 
horizontal curve of a minimum radius. Method 5 is perfect even when all vehicles are running at different speeds. Method 5is the 
preferred method suggested by AASHTO for superelevation distribution and side friction factor for all curves within radii larger than 
the minimum radius on rural highways and high speed urban roads. Method 5 has a curvilinear distribution involving an 
unsymmetrical parabolic curve (traditional curve) for the f distribution, which is the tangent of two legs that define Method 4. By 


subtracting the f rate from the design rate of (e + f ) for Eq. (1), the final superelevation distribution is achieved. 


3. SAFETY MARGIN 

Concept 

Nicholson (1998) defined the safety margin as the difference between the design speed (#) and the limiting speed (WL). The WL is 
the speed where f equals f,5, and the ¢ variable is the design superelevation. For a circular curve with radius R and design 
superelevation eg, the limiting speed is obtained by: 


‘ 


where FL= limiting speed (Km/h), andR= radius of a curve (m). Evidently, the design speed is less than or equal to the limiting speed 
(limiting speed equals the design speed if the design’sf equals f.,). Given that the value for fips, within the design speed tends to 
be less than the substantial available side friction, (limiting speed — design speed) presents a lower bound for the design safety 
margin. 

The difference in the safety margin along a highway segment with many horizontal curves is mainly due to the difference ine. It 
is tempting to distribute @ along a highway section where the difference for the safety margin can be small (improved design 
consistency) and the mean for the safety margin is greater (improved safety). Nicholson (1998) evaluated the safety margin using the 
AASHTO methods for the distribution off and e. 

Fig. 3 shows the mean and standard deviation of safety margin for all methods and optimization models. Method 2 had the 
smallest safety margin, while the largest safety margin was obtained with Method 3. Both Methods 2 and 3 did not fulfill the criteria 
for safety margin. The mean and standard deviation of safety margin for Methods 1 and 5 were between the values for Methods 2 
and 3. It was concluded that Method 1 (defined as the linear distribution) is better than Method 5 (defined as the curvilinear 
distribution). 

However, the outcome of Method 5 was better, dependent upon distributing @ (improperly) and utilizing a symmetrical parabolic 
curve (see Table 4). Additional analysis for Method 5 using an unsymmetrical parabolic curve defined by AASHTO showed that 
Method 5 worked best for high speed highway (Easa, 1999). 

Since the mean and standard deviation of the safety margin differ within the shape curve of the superelevation distribution, it is 
likely to locate the distribution to minimize the standard deviation of the safety margin, concerning to a limitation based on the 
minimum suitable safety margin. 
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IMPROVED SUPERELEVATION DISTRIBUTION FOR A SYSTEM OF CURVES 
To follow Easa’s (2003) approach, one must assume that the highway section under consideration has W circular horizontal curves 


with different radii, where the curves are divided into K groups with similar radii. Group = contains curves with radius R; , curve 
frequency g;, design speed ¥;, maximum side friction f mmx, design superelevation @;, design side friction f; and minimum radius 
R min; where t = 1,2,..... Both the minimum and maximum superelevations for the highway segment are designated with @,:, and 
fmagx -|n order to solve both thee; and ff, variables,R; , gz, Vy, F MAX; Smin, Sygx aNd R min; need to be known. The previous rotation 
symbolizes the wide scale case and that the maximum side friction and design speed (along with the minimum radius) differ within 
the curve radius of a certain speed environment, based on the Australian design standards (Rural, 1993). Using AASHTO standards, 
Fimax Renin, aNd V are similar for a certain highway section, and f max;= finge, Rming= Ri, and ¥;= ¥ for every i. 


Distribution using fixed curves 

The New EAU Method 

The new method is a proposed mathematical procedure for superelevation distribution based on AASHTO’s Method 5 and the EAU 
Curve Eq. (see Fig.2).This section describes the derivation of the proposed method. Side friction factors and superelevation 
distribution rate are all very important when designing appropriate horizontal highway alignments. Based on the laws of mechanics, 
the superelevation rate(@) needed for driving along a horizontal curve contributes to the following equation: 


ve 


0.01e = = 


=i (19) 


where, g = gravity constant (9.81 m/s), R = turning radius, w = vehicle's speed (km/h), and f = side friction factor. Practical design 
values exist for the upper limits of ¢ andg...in conjunction withfinc,e, taking into account many conditions associated with driving 
comfort, safety, pavement, traffic and weather. 

According to AASHTO, when traveling at the precise design speed the minimum turning radius (Ry _,) can be resolved such that 
both @yngxAN foc are chosen and shown in Eq.(2). This value functions as a maximum value for limiting both the side friction factors 


and superelevation rates from what is considered realistic for either comfort or operation by drivers. Therefore, utilizing radii larger 
than Rint, permits both superelevation rates and side friction factors design values that are below their upper limits. The curve 


parameters (Rap App §y.z»Ly.Lz) would be calculated according to AASHTO as follows: 
The radius at the point of intersection (Rzy) can be obtained from Eq. (20). 


~ LFF Gare 


where, PI is the Point of Intersection, with legs (1) and (2) for the f distribution parabolic curve, and Wy = running speed (km/h). The 
PI counterbalance from 1/R axis is shown in Fig. 2, and can be obtained from Eq. (21). 


(O01 BT 
Dapy = (SS SeeB) | 0.01 ence (21) 
The slope for leg (1),™, can be obtained from Eq. (22) 


G1 = Apr * Rp, (22) 


The slope for leg (2), gz, can be obtained from Eq. (23). 


Ear — 22 
s,- == (23) 


A is the algebraic difference betweengz and g,(A = gz-gy). The length of the first arc (Ly) and the second arc (Lz) and the total 
length of the curve (L}can be obtained by Eq. (24) to (26). 


Ly = 5 (24) 
i 4 

a (25) 

L=-L,+42, (26) 
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The curve parameter of the proposed EAU vertical curve equation is x and is given by: 


X= (27) 


1 
x 


The parameter that describes the unsymmetrical curve (Rzgy) can be obtained from Eq. (28): 


Reay = = (28) 


Rey is characterized as the ratio of the length of the shorter tangent (or the shorter arc when dealing with a traditional curve) to the 


total length of the curve. The rate of change of grades for the first arc (®,) and the second arc (rz) can be determined from Eq. (29) 


and (30), respectively. 


7, = L,<k, (29) 
ees (30) 


The value of is measured from the vertical axis. When x < 
determined as: 

a 
fa = 9.x + = (31) 


Note that ygyc= 0 for fy; where, f, = f distribution at any given point x < L/2; and 


4- oh (32) 


L/2, the first arc’s elevation for the EAU Curve; where, yy = ff can be 


whereg, = @ distribution at any given point x < L/2.When x>L/2, the second arc elevations for the EAU curve; where,yz = fcan be 


obtained as follows: 
fz = fmax - 92(k —x) + = — a = (33) 


Note that Yeyc= fimax for f,where fF, isthe f distribution at any given point x>L/2; and 


2" iaR fr (34) 


wheree ise distribution at any given point x> L/2. 


025 4 
> * (e+f)max 


@ fmax 
& fmax 


@ e/100max 


015 + 


€/100 mas 


¢€/100 or f or e/100 +f 


e/100 distribution 


0.0005 0.001 


j eal sn jo 


Figure 2 Development Procedures for Superelevation Distribution of the EAU Method 
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The New SAU Method 

The new SAU method is a proposed 
mathematical procedure for 
superelevation distribution based on 
AASHTO’s Method 5 and the SAU 
Curve Equation. This section describes 
the derivation of the proposed 
method. When designing a suitable 
horizontal alignment for highways, 
superelevation rates distribution (or 
the equivalent turning 
radii/curvatures) and_ side friction 
factors are very important. The laws of 
mechanics and superelevation rates 
influence the driver's ability to cope 
with turning along a horizontal curve. 
The SAU Method of superelevation 
determines the curve parameters such 
aS Rin Rep, Mpy, Ox, G2, 2y, 0.0, 
(A= 92-9), and x(Eq. 2, 20 to 27). 
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Using the SAU vertical curve equation the rate of change of slope (?gy_) can be obtained as: 


Move = (=) (Ly - 2 Ly) (35) 
t, which is constant, can be obtained as: 
t= )a- ba) (36) 


The value of x is measured from the vertical axis. When x < 1/Rgy, or *> 1/Rz, the arc elevations for the SAU Curve, as y = f, can be 
determined by: 


= Yove,2 , Fa 
f = gx + Wx? + ox (37) 


Note that Agye= 0 for F,where Raye =the elevation of PVC, f is the distribution at any given point x < 1/Ryy, or x> 1/Rp,, and 


v= 
aa f 


@= 


(38) 

where @ is the distribution at any given point x < 1/Rz,, or x> 1/Rp;. 

Safety Margin Computation 

WL is the speed limit at which f equals £,,,, and thee variable is the design superelevation. In a circular curve with radius R and the 


design superelevation e, the limiting speed from Eq. (18), and the constant is based on the constant development by AASHTO 
(1/127 = 0.00787). The safety margin can be calculated as: 


mm; = VI; - VD; Wi (39) 


whererm; is the safety margin for the group of curves; WD; is the design speed for the group of curvesi ; WL is the limiting speed; Sg; 
= NW; WL; is given by: 


The mean, variance, standard deviation (SD) and coefficient of variation(C¥) of the safety margin for the EAU Method can be 
calculated with the following equations: 


Mean = 4 —y (41) 
jen EE 
: _ ER, army iean)? 
Variance = “ya (42) 
SD= VVariance (43) 
CV = SD { Mean (44) 


The €¥, thus, has no units and is expressed as a decimal or a percentage. 

The coefficient of variation, CV, reflects safety and consistency. Therefore, a smaller CVindicates smaller SD (improved consistency) 
and larger mean (improved safety). The EAU method improves the design consistency by focusing on the system for horizontal 
curves in total, including curve radii and the frequencies of occurrence. The EAU method can be solved by mathematical software, 
such as Lingo. 

The SAU method is helpful when designing superelevation within the SAU curve equations, particularly if the highway segment 
contains many curves. To effectively improve highway consistency in the SAU method and to consider the SAU curve equations, it is 
important to first determine the cubic polynomial function and then the best superelevation distribution. This approach is extended 
for all horizontal curves on the highway. 

This information allows for representation of the safety margin along with horizontal curves in total which improves the design Or 
consistency. The horizontal curves on a highway segment are identified by 1, 2...., NW, where the value of WV is defined by the number S 
of curves. Then, the limiting speed (WL;) can be determined from Eq. (40), and the safety margin (%m;) can be determined from Eq. & 
(39). The mean, variance, standard deviation, and coefficient of variation of the safety margin can be calculated using Eq. (41) to (44). 
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The proposed methods improve the design consistency based on the safety margin, and can be solved by mathematical software, 
such as Lingo. 


New Distribution Using Parametric Cubic Model 
The general cubic curve designed for the f distribution along with a component of the aggregate analysis (Quadratic Model), which 


is provided by Easa (2003), forms the Parametric Cubic Optimization Model of the safety margin. Similar to Easa’s research (2003) 
regarding highway design consistency, it is attractive to reduce the standard deviation for the safety margin. Objective function 
includes variance as one of its nonlinear variables desirable. As an alternative, minimizing the maximum safety margin could 
minimize the variance. Thus, the objective function can be displayed as: 


Minimize Z = M (45) 


where Z=objective function along with M = maximum safety margin. The mean safety margin mean is equal to or greater than the 
desirable minimum value, 5,7, , stated as, 


mean > Sin (46) 
The other constraints related within the objective function is: 
my; < M, Wi (47) 


where m= safety margin for Curve Group; The mz; was utilized in in Eq.(39); ® = design speed for Curve Group #; along with the sum 
for g; = N. The constraints for Equation (39) include the limiting speed for Curve Group, WL; , as shown in Eq. (40). Design side 
friction factor f; and superelevation @;, can be calculated as: 


_ R a an a(® 7 af—enine) 
fi =f. Max; =} Amity ow lew Amin, (= (48) 
GOOTE? V 
8; = seu th '#;(49) 


The last constraints are related within the maximum side friction 7 mmax;along with the maximum and minimum superelevation 2» 


ANd Sree 

fi < f max;, ¥i (50) 
82 Opin, Vi (51) 
8:5 Omsy, Vi (52) 


The Mean, Variance, Standard Deviation, and Coefficient of Variation of safety margin can be calculated in Eq.(41)-(44). The objective 
function for Eq. (45) along with the constraints for Eq. (46) - (52) symbolizes a nonlinear optimization model. The optimization model 
summary is given in Table 2.The model contains (92) constraints, (6) decision variables, categorized by four variables having a 
subscript i (fj , @;¥L;, and mj) as well as two variables, ¢ and d. A summary of the nonlinear optimization model in Table 2 can be 
solved utilizing existing optimization software like LINDO (Schrage, 1991a; 1991b). This software deals with large scale problems, 
connecting 800 constraints and 4,500 variables. 


4. APPLICATION AND RESULTS 

When designing suitable horizontal alignment for highways, both the distribution of side friction factors and superelevation rates (or 
the equivalent turning radii/curvatures) are critical. The purpose of this section is to demonstrate the application of the methods 
using numerical examples. When evaluating a system of curves, all methods use the mean (x), standard deviation (SB) and 
coefficient of variation (CW) of the safety margin to determine the best superelevation distribution design. These variables can be 
utilized to measure both the safety margin and consistency. 


Data Selection 

The horizontal curve radius (centerline of the roadway) was used as the curve radius. The step by step calculations for determining 
the distribution of superelevation are presented in the following examples. Examples 1, 2, 3, and 4 are similar to those used by 
Nicholson (1998) and Easa (2003), allowing easy comparison of their methods with the proposed methods. The examples contain a 
single carriageway rural highway within rolling terrain that is expected to contain 20 circular horizontal curves with radii RB; = 100, 


www.discoveryjournals.org OPEN ACCESS 


Page O 


ARTICLE 


150, 200, 250, 300, 350 m and corresponding frequencies ofg; = 3, 4, 6, 4, 2, 1. The design guide of Australia recommends that a 90 
km/h speed environment have design speeds of¥;= 71, 76, 78, 80, 81, and 82 km/h. Then the @y.i, ANd Gyc~ Values are 0.02 and 0.10, 
respectively. The # mzvalues are 0.30, 0.28, 0.27, 0.26, 0.25, and 0.24. The minimum radius is determined from Eq.(2).The summary 


of the input data can be seen in Table 1. 


Table 1 
Input Data for Examples 1, 2, 3, and 4 


Curve RK @ J IR Cmax Smax Rein v? SR ein Rine Re Reo 
Ne. (m) (km/h) (km/h) 
l 100 3 71 63.90 0.1 03 99 0.40 0.02 124.04 132.31 396.93 
2 150 4 76 66.50 0.1 0.28 120 030 0.02 151.60 162.43 454.80 
3 200 6 78 68.25 0.1 0.27 129 0.24 0.02 165.19 177.43 479.06 
4 250 4 80 70.00 O11 0.26 140 0.20 0.02 179.98 193.82 503.94 
5 300 2 81 70.88 O.1 O25 148 0.17 0.02 191.34 206.65 516.61 
6 380 «(1 82 71.75 O.1 0.24 156 0.15 0.02 203.63 220.60 $29.45 


V=VD 


* Minimum mean of safety margin, Sin = 10 km/h 


For simplicity, the curve and path radius are considered equal. Assuming the minimum superelevation is 0.02, then the lowest radius 
utilizing minimum superelevation is determined using Eq. (13), where the lowest radii that depend solely on friction, Ry, , along with 


solely superelevation, R;,, is analyzed using both Eq. (7) and Eq. (10), respectively, providing the values highlighted in Table 1. The 
radii greater than 100m surpass the minimum radius; consequently there is flexibility in the proportions of the centripetal force given 


by side friction and superelevation. Data from Table 1 are used in the following sections to illustrate the results for the AASHTO, 
EAU, and SAU methods. 


Example 1 (AASHTO Methods) 

Table 2 illustrates the results for AASHTO’s five methods. The values for ¢ and f can be determined using Eq. (3) in Method 1; Eq. (5) 
and (6) in Method 2; Eq. (8) and (9) in Method 3; Eq. (11) and (12) in modified Method 2; and Eq. (14) to (17)in Method 5. For each 
radius, the limiting speed WL can be determined using Eq. (18) and the design safety margin WL —¥ can be anticipated. Utilizing the 
frequency information for each curve radius (Table 2)the mean, standard deviation, and coefficient of variation of the safety margin 
can be determined in each method(Table 2). 

Using Method 5, e, f, WL and WL-¥ values are 0.100, 0.297, 71.292 and 0.292(respectively) for curves with 100 m radius and 
0.078, 0.073, 117.257 and 36.951for curves with 350 radius (Table 2). The mean and standard deviation of the safety margin are 17.14 
and 10.71, respectively (Table 3). For the AASHTO methods the largest mean of the safety margin occurs when Method 3 or 5 are 
used (Table 3). However, the smallest standard deviation occurs when Method 2 is used. Therefore, none of the methods can be 
considered the best. Clearly, Method 5 provides a larger mean for safety margin (17.14 km/h) than Method 1 (13.8 km/h), Method 2 
(6.35 km/h) and modified Method 2 (8.54 km/h). The mean for safety margin in Method 5 (17.14 km/h) is only slightly smaller than 
that of Method 3 (18.44 km/h). Furthermore, the standard deviation of safety margin in Method 5 (10.71) is smaller than that of 
Method 3 (11.81) and larger than that of Method 1 (8.82), Method 2 (6.54), and modified Method 2 (8.01). Method 5 has the lowest 
coefficient of variance (0.63). The coefficient of variance was 0.64 for Method 1, 1.03 for Method 2, 0.64 for Method 3, and 0.94 in 
modified Method 2 (Table 3). Therefore, we conclude that Method 5 provides a good balance between the three measurements of 
safety and consistency and thus is considered the best method. 


Example 2 (EAU Method) 

Using the new EAU Method, e, f, WZ and WL—-¥ values are 0.100, 0.297, 72.031 and 1.031(respectively) for curves with 100 m radius 
and 0.078, 0.073, 120.121 and 38.121 for curves with 350 m radius (Table 2). The mean and standard deviation of the safety margin 
are 18.27 and 10.89, (Table 3). Table 3 compares this method with AASHTO’s methods in terms of the mean, SD and C¥ of the safety 
margin. 

The mean of the EAU Method (18.27) is higher (better) than the mean for Method 1 (13.80), Methods 2 (6.35), modified Methods 2 
(8.54) and Methods 5 (17.14). However it is slightly smaller than the mean of Method 3 (18.44) (Table 3). The standard deviation of 
the EAU method (10.89) is slightly smaller (better) than SD of Method 3 (11.81); however it is higher than SD of Methods 1, 2, 2 
modified and 5 (8.82, 6.54, 8.01 and 10.71, respectively) (Table 3). The coefficient of variation of the EAU method (0.60) has a margin 
of improvement from Methods 1, 2, 3, 2 modified and 5 (Table 3). The performance of different methods was examined 
USINGE cxf max, aNd Wy, where these variables were the same for all curves. The EAU method is superior to AASHTO’s methods. 
Discrete, and Quadratic model belongs to distribution @ and fF. 
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Table 2 Safety Margins for Difference Methods 


Radius (m) 


Method Variable 
100 150 200 250 300 350 
1 E 0.099 0.080 0.065 0.056 0.049 0.044 
F 0.298 0,223 0.175 0.146 0.123 0.107 


VL (km/h) 71.219 «82.805 92.227 100.185 106.792 112.477 
VL-V (knvh) 0.219 6.805 14.227 20.185 25.792 30,477 
2 E 0.097 0.023 0.000 0.000 0.000 0.000 
F 0.300 0.280 0.239 0.201 0.172 0.151 
VL(knv/h) 71.000 76.000 82.834 90.880 97.621 103.312 
VL-V (kivh) _0.000__0.000_4.834_10.880__16.621_21.312 
3 E 0.100 0.100 0.100 0.100 0.100 0.100 
F 0.297 0.203 0.139 0.101 0.072 0.051 
VL(knv/h) 71.292 85.104 96.968 106.939 115.507 122.966 
VL-V (knvh) 0.292 9.104 18.968 26.938 34.507 40.966 
2 modified E 0.097 0.023 0.020 0.020 0.020 0.020 
F 0.300 0.280 0.219 0.181 0.152 0.131 
VL(knv/h) 71.000 76,000 85.847 94.311 101.451 107.531 
VL-V (knvh) _ 0.000 0.000 7.847 14.311 20.451 25.531 
£ 0.100 0.097 0.091 0.087 0.082 0.078 
F 0.297 0.206 0.148 0.115 0.090 0.073 
VL(knv/h) 71.292 84.644 95.322 104.019 111.163 117.257 
VL-V (knvh) 0.292 8.792 17.843 24.926 31.484 36.951 
EAU E 0.100 0,098 0,093 0.088 0.083 0.078 
F 0.297 0.205 0.147 0.114 0.089 0.073 
VL(knv/h) 72.031 85.806 97.059 106.250 113.814 = 120.121 
VL-V (knvh) 1.031 9.806 19.059 26.250 32.814 38.121 
SAU E 0.100 0.100 0.096 0.089 0.082 0.076 
F 0.298 0.204 0.145 0.113 0.091 0.076 
VL(knv/h) 71.292 85,095 96.381 105.285 112.519 =—-118.526 
VL-V (knvh) 0.292 9.095 18.381 25.285 31.519 36.526 
Cubic Model E 0.100 0.068 0.044 0.032 0.024 0.020 
F 0.297 0.235 0.195 0.170 0.148 0.131 
VL(knv/h) 72.031 81.944 89.611 96.413 102.402. 107.617 
VL-V (knvh) 1.031 5.944 11.611 16.413 21.402 25.617 


a 


Table 3 
Comparison of Mean, Standard Deviation and Coefficient of Variation of Safety Margin of all Methods 


Safety margin (km/h) 


Coefficient of 
Element Mean Standard Deviation Variation 
__ AASHTO methods _ 
Method 1 13.80 8,82 0.64 
Method 2 6.35 6.54 1.03 
Method 3 18.44 11.81 0.64 
Method 2n? 8.54 8.01 0.94 
Method 5 17.14 10,71 0.63 
Fixed curve 
EAU Method 18.27 10.89 0.60 
Cubic *¢ 11.53 7.05 0.61 


"Based on Easa (2003). 
" The modified Method 2m is the same as Method 2 except that e¢,,,,. = 
0.02 instead of zero, 


© For the quadratic distribution, c=-0.000073. 
* For the cubic distribution, c=-0.004658587, 
" For the cubic distribution, c=-0.0002452192. 


Example 3 (SAU Method) 

Using the SAU method, the ge, f, VL and VL-¥ values are 0.100, 0.298, 71.292 and 0.292 (respectively) for curves with 100 m radius 
and 0.076, 0.076, 118.526 and 36.526 for curves with 350 m radius (Table 2). The mean, standard deviation and coefficient of 
variation of the safety margin are 18.39, 10.88 and 0.59, respectively (Table 3). Table 3 illustrates that SAU is a better method than 
the AASHTO methods in terms of the mean, 5D and €¥ of safety margin. The mean of the SAU method has a margin of 
improvement from Methods 1, 2, 2 modified and 5; however, it is slightly smaller than the mean in Method 3 (Table 3). The standard 
deviation of the SAU method has a margin of improvement from Method 3; however it is slightly bigger than that of Methods 1, 2, 2 
modified and 5 (Table 3). The coefficient of variation of the SAU Method has a margin of improvement from Methods 1, 2, 3, 2 
modified and 5. 
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The superelevations differ drastically between the improved method and AASHTO’s Method 5. While AASHTO’s superelevation 
distribution utilizes a parabolic formula for an unsymmetrical curve, the SAU method uses a cubic formula for the SAU curve. 

The performance of different methods was examined using the same@mox, fing, aNd Wp variables for all curves. The SAU method 
is superior to the AASHTO method, EAU method, Discrete, and Quadratic models in regards to distribution @ and f. 


Example 4 Parametric Cubic Model (Cubic distribution of fF) 
The results of the optimization model are given in Table 2. Using the LINGO software based on the Global optimal solution, the 
objective function was obtained to be 25.617, and the parameters of the optimization model determined that ¢ = —0.004658, and d 


= 0.000245. The range of the safety margin is from 1.03 to 25.62 km/h (Table 2) within an 11.53 km/h mean, 7.05 Standard Deviation 
and 0.61 Coefficient of Variation (Table 3). The optimization model presented here is the improved solution based on the CV of the 


safety margin for this new model is better than that of AASHTO methods. A comparison for both the superelevation @; and the 
safety margin mj; for AASHTO Method 5, along with the optimization model is illustrated within Table 3. The superelevations differ 
drastically between the optimization model and AASHTO Method 5. While AASHTO superelevation distribution utilizes a parabolic 
formula for an unsymmetrical curve, the optimization model utilizes a cubic formula of f distribution. The optimization model 


minimizes the variation in safety margin of a specific highway segment. 
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Figure 3 Mean and Standard Deviation of Safety Margin for All Methods and Optimization Models 
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Figure 4 Coefficient of Variation of Safety Margin for All Methods and Optimization Models 


5. DISCUSSION 

This paper presents two proposed methods and an optimization model, the EAU, SAU Methods and Parametric Cubic Model, 
respectively. A comparison of superelevation rate ;, side friction factor f, speed limit WL; and the safety margin m; of AASHTO 
Methods and the optimization model are shown in Table 2. For the counterbalance of the mean, standard deviation and the 
coefficient of variation of the safety margin, the results demonstrated here appear in relation to the curvilinear distribution 
belonging to Method 5 that is recommended by AASHTO. Table 3 shows all the AASHTO, the EAU, SAU Methods along with the 


www.discoveryjournals.org OPEN ACCESS 


Paseo 


ARTICLE 


Discrete, Quadratic and Parametric Cubic Optimization Models. The table shows that the SAU Method improves the ¢ and f 
distribution, factors that depend on the frequency and radii of the curves, roadway restriction, safety margin and design consistency. 
Fig. 3 shows the mean and SD of safety margin for all methods and optimization models. Comparing the EAU with all other 
methods and models, it observes a counterbalance for the mean of the safety margin where the EAU method is superior to 
AASHTO's methods, and the Discrete, Quadratic and Parametric Cubic Models (a greater mean translates to better safety). However, 
EAU is inferior to the SAU method in improving design consistency based on the safety margin. Considering the Parametric Cubic 
Model, it is the best in order of ¢ and f distribution (Table 2) and it is superior to AASHTO Methods and the Quadratic model; 
however, it is slightly inferior to the SAU and EAU Methods and Discrete Model in maximizing design consistency based on the 
safety margin. Additionally, based on the coefficient of variation values (Fig. 4), EAU and SAU methods are superior to all AASHTO 
methods and optimization models (a lower coefficient of variation translates to higher design consistency). 


6. CONCLUSION 

This paper presented an analysis of superelevation design using two new methods. The proposed methods and the optimization 
model for superelevation design are easy to apply and provide superelevation rates that are comparable with those of AASHTO’s 
Method 5. As explained earlier, Method 5 signifies a mathematical analysis which has little consideration for speed variation, as well 
as entailing an extensive process to attain the superelevation distribution. For the AASHTO methods, inconsistency in safety margin 
could have been due to utilizing curves with radii greater than the minimum radius. In some situations, the variation in safety margin 
is due to superelevation and side friction factor. The system of curves based on the EAU and SAU methods offers a reasonably high 
mean of safety margin (i.e., higher than Method 5) and a low coefficient of variation. The aforementioned translates into a higher 
consistency level for the safety margin (i.e., higher than Method 5). This paper presents an optimization model which considers the 
mean, standard deviation, along with the coefficient for variation for the probability distribution for each variable (@ and f) that are 
different from the extreme values. The optimization model also accounts for the connection between variables. This paper proposes 
a new technique to improve highway design consistency following the principles of safety margin. The safety margin is defined by 
the difference between the design speed and the maximum limiting speed. The advantages of the proposed techniques are as 
follows: 


1. They eliminate the requirement for trial-and-error in evaluating all AASHTO methods to find the superelevation design for any 
known set of data for highway curves. Moreover, these methods directly find the best superelevation taking into account the 
complete design space of the curve (of which AASHTO curves are a division). 

2. Both the EAU and SAU methods are valuable when the highway section has numerous horizontal curves. The aforementioned 
increases ad hoc expectancy for drivers on a specific highway. This study is theoretically similar to AASHTO Method 5, which 
offers an even distribution of # for highway curves, and succeeds analysis by Nicholson (1998) and Easa (1999; 2003). 
Consequently, this analysis presents a clear method for deliberation of the safety margin when designing superelevation to 
improve highway consistency. 

3. Since the methods being used involve unsymmetrical vertical curve equations, the most compatible solution is assured. This 
method can also take care of long highway sections and can effortlessly solve issues relating to any useful number of highway 
horizontal curves using available powerful mathematical optimization software. 

4. In consequence, the optimization model determines the superelevation rate and side friction factor that minimize the objective 
function of a certain highway segment. 

5. Highway design consistency based on the safety margin approach assists that support greater consistency in selecting the 
design speed. The method proposed in this paper will definitely assist the safety margin approach. 


The proposed methods provide a superelevation design that allows a larger safety margin for most drivers. The information 
provided in this paper is helpful for both highway evaluators and designers, especially when considering design sufficiency. 
Integrating the safety index within the turning radius design, and the superelevation distribution promises a thorough proposed 
methods analysis that are built within the design methodology. This study considers the superelevation distribution of horizontal 
curves within highway design in relation to the effects of several parameters; based on the numerical examples provided, it was 
determined that the new EAU, SAU Methods and the Parametric Cubic Model were superior to AASHTO Methods, due to 
improvements in highway design consistency. Future work needs to consider the reliability index as a possible improvement to 
safety margin. 
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Appendix : Parametric Cubic Model derivation (Cubic distribution of f) 
Consider the general cubic curve designed for the f distribution, exposed by dash line within Fig. 1(b). The Eq. for the curve is 
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f=a+b(=) e(2) +a(2) (53) 

At R= 00, f=0, and consequently @=0. Furthermore, at R=Rynin » F= Frmax: 

fos O() el) dG) 8 

Based on Eq. (54) 

b = fmax®min -¢(;-—-)- 4 (4) (55) 

Substitute b from Eq. (55) in Eq. (53) 

f- = —c(=)-a (2)1@+ c(2) +a ey (56) 
Eq. (56) can be re-written as follows: 

pe Measiase «(Sam () -2 (“Ee )G) sn 


By replacing the parameter (1/R) with (1000/R) in order to keep away from the scaling problem to get the optimum solution for 
Curve Group i (note that a=0) 


fictmen Se CRN EEC) 6 


where ¢ and @ = indefinite variables (decision variables). The decision variables ¢ and @ might be negative or positive. For d=0 


Eq.(58) is quadratic. 


Notation: 
The following symbols are used in this paper: 


A = algebraic difference in grad; 

a,b,c. = parameters of cubic function of f distribution; 
C¥ = coefficient of variation of the safety margin; 

@, f = superelevation and side friction factor, respectively; 
@; = design superelevation of Group i; 

Smin = Minimum superelevation; 

fmqx = Maximum permissible design superelevation; 
@z, = superelevation utilizing Method 2 (modified); 
fimax = Side-friction coefficient; 

fi = design side friction of Group i; 

fmax; = maximum side friction of Group i; 

g = acceleration of gravity; 

9, = slope for leg 1; 

92 = slope for leg 2; 

i = index for curve group, where i=1,2....,K; 

j = index for set of curves, f=1,2,...,N; 

Rtg; = Pl counterbalance from the 1/R axis; 

Rgye = elevation of PVC for SAU curve; 

K = number of curve groups; 

L, = length for the first arc; 

L, = length for the second arc; 

L = total length for the curve; 
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M = maximum safety margin; 


mj; = safety margin of group i; 

Mean = mean of the safety margin; 

NN = total number of curves in highway segment; 

g; = curve frequency for Group i; 

R = curve radius; 

R; = curve radius for Group i; 

Roig, = Minimum curve radius; 

Rz, = smallest potential radius relying by superelevation only; 
Ry = smallest potential radius relying by side friction only; 
Riune= smallest potential radius utilizing minimum superelevation; 
Ry, = radius of at the point of intersection, Pl; 

Regy = the ratio for the length of the shorter tangent for EAU curve; 
¥, = rate of change of grad for the first arc for EAU curve; 

¥z = rate of change of grad for the second arcfor EAU curve; 
pyc = rate of change of slop for SAU curve; 

Swi, = Minimum allowable mean safety margin; 

5D = standard deviation of the safety margin; 

t = the third derivative of y for SAU curve; 

W = design speed; 

¥; = design speed for Group i; 

WL = limiting speed; 

Vg= running speed, km/h; 

Variance = variance of the safety margin; 


Yave = elevation at BVC point for EAU curve;Vvzey_e = elevation at EVC point for EAU curve; and 


Z = objective function. 
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